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(57) Abstract: The invention concerns a method for determining a bone mechanical strength from a digitized two-dimensional 
image, obtained by imaging, characterized in that it consists in correlating the bone mineral density determined from said two- 
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a partir de cette image 2D par tout mo yen adapte a ce type d' image et un parametre de structure obtenu a partir de la meme image 2D. 
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PROCEDE DE DETERMINATION D'UNE STRUCTURE EN 3D A PARTIR 
D'UNE IMAGE EN 2D, NOTAMMENT LA STRUCTURE D'UN OS 

La presente invention concerne un procede de determination d'une structure en 
3D a partir d'une image en 2D, notamment celle d'un os. 

L'osteoporose est une maladie qui affecte en grand nombre surtout les femmes 
apres la menopause vers Tage de 50 ans sachant que cette maladie peut affecter 
5 toute personne et a tout age. Elle est caracterisee par une faible masse osseuse 
et une deterioration du tissu osseux. Ces degradations conduisent a des risques 
importants de fractures particulierement de la hanche, de la colonne vertebrate 
et du poignet. 

II existe bien sur un certain nombre de facteurs de risques mais ceci ne suff it 
10 pas pour determiner la probability pour une personne de se f racturer un element 

de son squelette sachant qu'il n'y a aucun symptome prealable. 

Cette maladie est courante dans la population et devrait aller en s'etendant du 

fait de la sedentarisation notamment et du vieillissement de la population. 

Or le traitement des personnes aff ectees apres fracture est tres long et tres 
15 couteux car ce sont des soins a long terme. Les consequences sont souvent 

graves car conduisant a des invalidites voire a la mort des patients. 

Aussi est-il necessaire de prevenir ce type d'affection mais d condition de 

determiner le risque de fracture, ce qui est du seul ressort du praticien. 

Si Ton arrive a predire les risques, le praticien dispose alors de traitements par 
20 hormonotherapie et par prise d'autres produits de la pharmacopee comme la 

calcitonine ou les biphosphonates, en plus de conseils d' hygiene de vie pour 

tenter d'enrayer cette affection. 
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II existe des appareillages pour determiner la densite osseuse appeles 
densitometres. 

On peut penser que des lors que Ton conncut la densite osseuse, en comparant 
avec des courbes etalon, le praticien peut alors determiner les risques et etablir 
5 un diagnostic adapte et prescrire un traitement satisf aisant. 
C'est ce qui est couramment utilise. 

Par la demande de brevet US 5 774 520, on sait que la probability d'une 
fracture est directement liee a la densite osseuse. 

Les densitometres disponibles dans le commerce utilisent les rayons X ou les 
10 rayons gamma de faible energie. L'absorption etant proportionnelle a la densite 
de I'os, on peut done penser a correler les deux pour avoir un resultat 
satisfaisant permettant des comparaisons. Ce brevet prevoit un travail avec 
deux sources d' emission de photons en sorte de traiter la dispersion engendree 
par le fait que les matieres osseuses peuvent etre considerees comme 
15 homogenes mais que les autres matieres liquides, les tissus musculai res, les 
cartilages ne peuvent etre consideres comme une matiere unique, homogene. 
Un autre brevet US 6 385 283 utilise la densite mais combine cette mesure avec 
des images permettant de determiner les risques de fracture. Ces images sont 
realisees sur la colonne vertebrate du patient et I'operateur determine la 
20 presence d' amorces de fractures. De plus, comme dans Tart anterieur connu, 
ces informations sont comparees avec des donnees acquises au fur et a mesure 
sur de tres nombreux cas. Ce diagnostic est aff ine par le praticien a partir des 
antecedents du patient et par introduction des facteurs de risques 
supplementaires. 

25 La demande de brevet WO 86/07531 propose de realiser une image d'un os 
donne, par exemple le calcaneum, a un age auquel le patient presente 
normalement sa plenitude de masse osseuse. Ensuite, le procede consiste a 
realiser des images successes de ce meme os pour permettre au praticien de 
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realiser des comparaisons et de determiner revolution afin de conclure sur le 
degre de risque de fracture. 

Toutes ces methodes et autres dispositifs sont bases sur la seule mesure de la 
densite de la matiere osseuse. 
5 Or pour que le praticien arrive a un diagnostic satisfaisant, il faut qu'il ait des 
informations autres que la seule densite. 

La definition de I'osteoporose est donnee par I' Organisation Mondiale de la 
Sante : "desordre caracterise par une faible masse osseuse et des alterations 
de la microarchitecture du tissu osseux, conduisant a augmenter la fragilite de 

1 0 I'os et en consequence a augmenter le risque de fracture " 

C'est ainsi qu'il existe des analyses de matiere osseuse in vivo plus sophistiquees 
a partir d'analyses de sang et d'urines faisant intervener des marqueurs osseux 
mais ces analyses sont couteuses, contraignantes pour le patient et, de toutes 
les fagons, donnent des indications de composition de la matiere mais pas 

15 d' indications sur sa structure. 

Or, on peut comprendre que la matiere osseuse presente une resistance 
differente suivant qu'il y a plus ou moins de matiere en quantite ce qui resulte 
de la mesure de densite. Par contre, pour une meme densite osseuse done pour 
une meme quantite de matiere, suivant I' architecture de la structure osseuse, la 

20 resistance mecanique peut varier grandement. 

Ce qui est utile pour un praticien, ce sont les proprietes mecaniques de I'os d'un 
patient induisant directement les risques de fracture. Quand il y a diminution de 
la resistance osseuse, il y a conjugaison de deux facteurs, la baisse de la densite 
minerale osseuse DMO et la desalteration de la structure osseuse. 

25 A ce jour, il n' existe pas de moyens permettant de determiner la structure 
osseuse. On peut utiliser de I' imager ie par resonance magnet ique ou le scanner 
mais ce sont des examens lourds et couteux, surtout s'il faut effectuer un suivi 
longitudinal des patients. 
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De plus, la capture d'images en 3D a partir de moyens d'imagerie est pour 
['instant limitee par les capacites de ces moyens car la resolution spatiale est 
superieure aux dimensions de la plupart des travees de liaison permettant de 
realiser la connectivity, c'est-a-dire les liaisons entre les noeuds du reseau en 

5 trois dimensions. 

Le procede selon la presente invention est maintenant decrit en detail, afin de 
permettre par I' analyse d'une image 2D d'estimer les parametres mecaniques. 
Selon I'invention, le procede de determination de la resistance mecanique d'un os 
selon T invention, a partir d'une image 2D numerisee, obtenue par imagerie, se 

10 caracterise en ce que Ton realise une correlation entre la densite minerale 
osseuse DMO determinee a partir de cette image 2D par tout moyen adapt e a ce 
type d'image et un parametre de structure obtenu a partir de la meme image 2D. 
Plus particulierement, on recourt a une correlation de type exponentielle. 
On utilise cette correlation associant la densite minerale osseuse et ledit 

15 parametre de structure pour determiner la contrainte ultime C u de I'os. 

De f agon plus precise, on determine le parametre de structure a , obtenu par la 
succession des etapes suivantes : 

a) choisir au hasard un pixel de I' image 2D qui est a un niveau de gris h(0), 

b) choisir une droite issue de ce point ayant une direction egalement determinee 
20 au hasard, 

c) se deplacer d'une distance r sur cette droite, h(r) etant le niveau de gris de 
ce nouveau point, 

d) calculer la variance des niveaux de gris avec la formule : V(r) = [h(r) - h(0)] 2 , 

e) tracer la courbe associee a V(r) en echelle log-log, et 

25 f ) determiner la pente de cette courbe log-log qui represente ledit parametre a. 
Afin d'ameliorer la precision, on repete les etapes a) a d) un nombre de fois 
suffisamment grand pour faire converger la fonction variance V(r), moyennee 
sur I 1 ensemble des iterations. 
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Selon une autre caracteristique de I' invention, on realise une correlation entre la 
densite minerale osseuse, DMO, obtenue a partir de cette image 2D et ledit 
parametre a evalue a partir de la meme image 2D suivant le modele 
mathematique : 

5 C u ' = b 0 + bi* exp (b 2 * DMO) * a 

- avec bo, bi, bz qui sont des coefficients obtenus par regression non 
lineaire et C u ' la prediction de la contrainte ultime C u de I'os. 
On determine une correlation entre le parametre a et un parametre 3D du 
reseau trabeculaire de I'os et un parametre 3D peut etre la densite de 
10 connectivity Xv 

Pour supporter cette description, des dessins sont annexes et les figures qu'ils 
presentent, illustrent la description et montrent essentiellement des resultats 
obtenus. 

- figure 1 : courbe des resultats avec un modele lineaire, C u en f onction de C u « 

15 - figure 2 : courbe des resultats avec un modele exponential, C u en f onction de 

Cu' 

- figure 3 : representation graphique de C u en f onction de Xv 

- figure 4 : representation de f s en f onction de Xv 

- figure 5 .courbe de la f onction V(r) 

20 - figure 6 : determination du parametre a a partir de la courbe log-log de la 
f onction V(r) 

- figure 7 : courbe DMO en f onction de f s 

- figure 8 : courbe de a en f onction de Xv 

- figure 9 : courbe DMO en f onction de a 

25 - figure 10 : modele lineaire C u en fonction de DMO 

- figure 11 : courbe de C u en fonction C u \ 

L 1 etude porte sur un os trabeculaire ou spongieux qui constitue environ 20% de 
la matiere osseuse. L'os cortical qui Tentoure assure la rigidite de Tos complet 



WO 2004/057512 PCT/FR2003/003768 

6 

et done du squelette tandis que I'os trabeculaire assure la reprise et la 
resistance aux forces de compression. 

Selon T invention, on considere que la structure trabeculaire comprend des 
travees qui sont interconnectees et la resistance mecanique de I' ensemble 
5 resulte de cette connectivity e'est-a-dire du nombre de boucles fermees. La 
resistance mecanique resulte aussi du degre de mineralisation de ces travees, 
mineralisation qui est un parametre important. 

Le procede consiste a analyser une image numerisee obtenue dans le cas present, 
a partir d'une emission de rayons X. Cette image est une projection sur une 
10 surface d'une structure 3D ayant ete traversee par les photons emis. 

Chaque pixel qui constitue 1 1 image peut etre analyse independamment mais alors 
on ne peut determiner que la densite minerale osseuse DMO si bien qu'il convient 
non seulement d' analyser un pixel de fagon isolee mais de les analyser les uns par 
rapport aux autres. 

15 Le spectre du rayonnement doit etre le plus monochromatique possible afin 
d'eviter une trop forte derive dans I'energie des photons. 

En effet, si I'energie varie la penetration variera et donnera des nuances de gris 
sur T image pouvant induire des erreurs en faisant croire a une absorption la ou il 
n'y en a pas. 

20 A partir de cette image obtenue avec un spectre etroit, on sait que chaque 
volume elementaire subit une action directe des photons, perpendiculairement 
mais aussi une action indirecte des rayons incidents diffuses par les volumes 
elementaires voisins apres qu'ils aient ete eux-memes tras/erses. 
Afin de supprimer la part de diffusion, on peut utiliser de preference deux 

25 energies distinctes. On dispose de 2 equations a 2 \r\conr\ues et on peut done 
deduire la part d'attenuation reelle qui interesse directement le present 
procede et la part resultant de la diffusion. 
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Si Ton recherche la contrainte ultime C u (MPa), il faut tenir compte des 
parametres suivants : 

- Vo*. volume d' analyse (mm 3 ) 

- V s : volume du solide (mm 3 ) 

5 - f s : fraction volumique de solide 
avec la relation f s = V s / Vo 

- Po : nombre d'amas de la partie solide (c'est-a-dire I' ensemble des parties 
connectees ou disconnectees), ce nombre est generalement egal a 1. 

- p2 nombre de surfaces internes, c'est-a-dire les trous resultant de la micro 
10 porosite interne des travees or a I'echelle de travail, cette porosite ne se voit 

pas, ce nombre est generalement egal a 0. 

- Nep nombre d'Euler-Poincarre 

- % : connectivity 

avec la relation x = Po + P2 - Nep 
15 - xv *• densite de connectivity (mm -3 ) 
avec la relation Xv = X / Vo 



Le procede selon la presente invention vise a permettre de correler deux 
parametres mesurables et de trouver la relation qui les lie. 
20 1 / Si Ton utilise un modele lineaire tel que : 

Cu = ao * fs + ai 

- C u ' etant une valeur predite de C u 

- ao et ai sont des coefficients de regression lineaire, 
il s'agit la d'une simple correlation de C u avec f s . 

25 La courbe obtenue est representee sur la figure 1 et Ton constate que la 
dispersion est tres importante. Si Ton se place par exemple a f s constant, 
I'ecart des valeurs de C u est important. 
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On ne peut retenir ce modele lineaire dont les resultats sont insuffisants en 
precision mais qui peut neanmoins servir de comparatif . 

2/ Si Ton utilise un modele exponentiel tel que : 
5 C u '= bo + bi* exp(b 2 * f s ) * Xv 

Avec bo, bi, bz qui sont des coefficients obtenus par regression non lineaire 
(methode des moindres carves generalisee) et C u ' la prediction de la contrainte 
ultime C u de I'os. 

La dispersion est fortement limitee, voir figure 2. C u ' est estime a partir de f s 

10 etxv. 

Ces variations sont representees sur la figure 3. 

Si f s et Xv augmentent, C u augmente. Ceci est parfaitement naturel, si la fraction 
volumique de solide augmente et la connectivity augmente, alors la contrainte 
ultime augmente. 

15 Si on fixe Xv alors C u est directement fonction de f S/ de fagon exponentiel le. A 
connectivity constante, la contrainte ultime augmente en fonction de 
I' augmentation de la densite osseuse. 

Si on fixe f s alors C u est directement fonction de Xv, de fagon lineaire mais la 
contrainte ultime C u diminue lorsque la connectivite augmente. 
20 Si la fraction solide est constante, on constate que la resistance ultime C u 
diminue lineairement en fonction de la connectivite. 

Plus on a de noeuds dans une structure et plus cette structure devient fragile, 
ceci de f agon surprenante, contrairement a un prejuge bien etabli. 
On peut alors tracer le modele de representation de la figure 4 qui montre les 
25 valeurs de C u . 

Pour la determination de la connectivite, il faut done trouver un parametre qui 
suive une meme loi et qui soit le plus independant possible de la densite osseuse 
pour ne pas etre influence. 
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On pourrait realiser un calcul a partir de deux parametres mesurables comme la 
fraction volumique osseuse et la connectivity, ces deux parametres etant 
obtenus notamment par Imagerie a Resonance Magnetique. 

Par centre, I'obtention des valeurs par cette methode met en oeuvre des 
5 appareillages couteux, peu disponibles, ce qui conduit a des predictions d'un prix 
de revient eleve interdisant un suivi longitudinal et done une diffusion massive 
aupres des patients. 

II faut done pouvoir trouver deux valeurs mesurables par exemple a partir d'une 
simple densitometrie, qui est, elle, tout a fait abordable et qui peut etre 
10 renouvelee regulierement. En effet, des densitometries sont deja realisees de 
fa^on massive pour de nombreuses applications. La diffusion aupres des patients 
peut etre tres importante. 

Le procede selon la presente invention consiste a determiner ces deux 
parametres mesurables comme la fraction volumique osseuse et la connectivity, 

15 ou plus exactement des parametres qui pourraient leur etre correles de fagon 
suf f isamment etroite pour donner des resultats satisfaisants. 
Le premier parametre rapidement mesurable et fiable est la densite minerale 
osseuse. Des appareils ont ete developpes pour effectuer ces mesures de fagon 
fiable et reproductible. Si Ton trace la courbe DMO en fonction de f s , on 

20 constate que la dispersion est faible et que Ton peut utiliser DMO en lieu et 
place de f s .. Voir figure 7. Le premier parametre est obtenu. 
Le second parametre est appele a et determine a partir d'une image 2D 
numerisee par exemple obtenue avec un rayonnement X. 

Pour determiner ce nouveau parametre a, on choisit tout d'abord au hasard un 
25 pixel de I' image 2D qui est a un niveau de gris h(0). 

On choisit ensuite une droite issue de ce point ayant une direction determinee 
egalement au hasard. 
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On se deplace d'une distance r sur cette droite, h(r) etant le niveau de gris de 
ce nouveau point. 

On calcule alors la variance des niveaux de gris, ce qui donne : 

V(r)=[h(r)-h(0)] 2 

5 Le processus est iteratif , avec une repetition de ces etapes un nombre de fois 
suff isamment grand, et aleatoire pour determiner a. On fait ainsi converger la 
fonction variance V(r) moyennee sur I' ensemble des iterations. Ceci permet 
d'obtenir par calcul un ensemble de points et la courbe associee est representee 
sur la figure 5. 

10 Ainsi, plus on determine de points ayant un meme niveau de gris a proximite du 
pixel choisi, plus on se trouve dans la matiere et done plus on est homogene et 
plus on doit trouver de noeuds de connexion. 

Si Ton trace la courbe log-log de cette fonction, figure 6, on trouve sur les cinq 
premiers points en V occurrence, un comportement lineaire et la pente de cette 
15 droite est le coefficient denomme a et propre a la presente invention. 

Si Ton trace maintenant la courbe a en fonction de la connectivity volumique Xv, 
figure 8, on constate qu'il y a une correlation suffisante pour que Ton puisse 
remplacer Xv par a. 

On peut deduire le diagramme de la figure 9 de ces courbes precedentes, 
20 diagramme qui montre la relation entre DMO, a et C u . 

DMO et a sont deux parametres directement mesurables et/ou pour le moins 
calculables a partir d'une meme image en 2D. 

On constate sur ce diagramme que lorsque la densite minerale osseuse reste 
constante, la contrainte ultime diminue lorsque a augmente. 
25 Ceci est bien equivalent au modele exponential precedemment indique avec la 
fraction osseuse en fonction de la connectivity, voir figure 4, sauf que la 
connectivity Xv n'etait pas mesurable directement alors que a est calcule a 
partir de I' image 2D, done quantifiable. 
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On peut done ecrire la meme relation mais avec les parametres determinates d 
parti r de la meme image : 

Cu = b 0 + bi* exp(b 2 * DMO) * a 
Cest la representation de la figure 9. La courbe de la figure 11 caracterise ce 
5 modele exponential de C u ' en fonction de a. On constate une faible dispersion 
comparee a celle du modele lineaire dont la courbe est mentionnee a la figure 10. 
Ceci est un exemple non limitatif qui montre la possibility de lier la densite 
osseuse dont la valeur est determinee a partir d'une image 2D, a un parametre 
de structure determine a partir de cette meme image 2D pour obtenir la 
10 contrainte ultime C u . 

On constate qu'il est possible de determiner une correlation entre le parametre 
a evalue a partir de V image 2D et une caracteristique 3D du reseau trabeculaire 
de I'os comme la densite de connectivity Xv. 

On peut ainsi donner un exemple obtenue a partir d' images osseuses femorales 
15 obtenue absorptiometrie a double energie qui remplace une image obtenue par 
rayons X ou Gamma de faible energie. 

Ces images sont celles Tune d'une jeune femme de 17 ans et I'autre d'une 
femme agee de 69 ans. Les images 2D ainsi obtenues sont analysees suivant le 
procede de T invention. 
20 On constate a la mesure de la DMO les memes valeurs obtenues directement par 
le systeme densitometrique, dans les deux cas : 0,730 g/cm 2 dans la zone 
trochanterielle. 

Les parametres a calcules sont respectivement 0,56 pour le sujet jeune et 0,71 
pour le sujet age, montrant ainsi une difference notable et le parametre a le 
25 plus faible correspond a la meilleure resistance mecanique de I'os C U/ celui du 
sujet le moins age. 
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REVENDICATTONS 

1. Procede de determination de la resistance mecanique d'un os a partir d'une 
image 2b numerisee, obtenue par imagerie, caracterise en ce que Ton realise une 
correlation entre la densite minerale osseuse DMO determinee a partir de cette 
image 2D par tout moyen adapte a ce type d'image et un parametre de structure 

5 obtenu a partir de la meme image 2D. 

2. Procede de determination de la resistance mecanique d'un os selon la 
revendication 1, caracterise en ce que Ton recourt a une correlation de type 
exponentielle. 

3. Procede de determination selon la revendication 1 ou 2, caracterise en ce que 
10 la correlation associant la densite minerale osseuse et ledit parametre de 

structure est utilisee pour determiner la contrainte ultime C u de I'os. 

4. Procede de determination de la resistance mecanique d'un os selon Tune 
quelconque des revendications precedentes, caracterise en ce que Ton determine 
le parametre de structure a obtenu par la succession des etapes suivantes : 

15 a) choisir au hasard un pixel de I' image 2D qui est a un niveau de gris h(0), 

b) choisir une droite issue de ce point ayant une direction egalement determinee 
au hasard, 

c) se deplacer d'une distance r sur cette droite, h(r) etant le niveau de gris de 
ce nouveau point, 

20 d) calculer la variance des niveaux de gris avec la formule : V(r) = [h(r) - h(0)] 2 , 

e) tracer la courbe associee a V(r) en echelle log-log, et 

f ) determiner la pente de cette courbe log-log qui represente ledit parametre a. 

5. Procede de determination de la resistance mecanique d'un os selon la 
revendication 4, caracterise en ce que Ton repete les etapes a) a d) un nombre 
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de fois suffisamment grand pour faire converger la fonction variance V(r), 
moyennee sur I' ensemble des iterations. 

6. Procede de determination de la resistance mecanique d'un os selon la 
revendication 4 ou 5, caracterise en ce que Ton realise une correlation entre la 

5 densite minerale osseuse, DMO, obtenue a parti r de cette image 2D et ledit 
parametre a evalue a partir de la meme image 2D suivant le modele 
mathematique : 

Cu = b 0 + bi* exp (b 2 * DMO) * a 
- avec bo, bi, bz qui sont des coefficients obtenus par regression non 
10 lineaire et C u ' la prediction de la contrainte ultime C u de Tos. 

7. Procede de determination de la resistance mecanique d'un os selon Tune 
quelconque des revendications 4, 5 ou 6, caracterise en ce que Ton determine 
une correlation entre le parametre a et un parametre 3D du reseau trabeculaire 
de Tos. 

15 8. Procede de determination de la resistance mecanique d*un os selon la 
revendication 7 caracterise en ce que le parametre 3D du reseau trabeculaire 
de I'os est la densite de connectivity Xv 
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